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ABSTRACT 

Interplanetary dust (IPD) scatters solar radiation which results in the zodiacal light that dominates 
the celestial diffuse brightness at optical and near-infrared wavelengths. Both asteroid collisions 
and cometary ejections produce the IPD, but the relative contribution from these two sources is 
still unknown. The Low Resolution Spectrometer (LRS) onboard the Cosmic Infrared Background 
ExpeRiment (CIBER) observed the astrophysical sky spectrum between 0.75 and 2.1 /im over a wide 
range of ecliptic latitude. The resulting zodiacal light spectrum is redder than the solar spectrum, 
and shows a broad absorption feature, previously unreported, at approximately 0.9 /im, suggesting the 
existence of silicates in the IPD material. The spectral shape of the zodiacal light is isotropic at all 
ecliptic latitudes within the measurement error. The zodiacal light spectrum, including the extended 
wavelength range to 2.5 /im using Infrared Telescope in Space (IRTS) data, is qualitatively similar to 
the reflectance of S-type asteroids. This result can be explained by the proximity of S-type asteroidal 
dust to Earth's orbit, and the relatively high albedo of asteroidal dust compared with cometary dust. 
Subject headings: infrared: diffuse background — interplanetary medium — zodiacal dust 



1. INTRODUCTION 

The astrophysical sky brightness is dominated by the 
zodiacal light which comprises scattered sunlight by in- 
terplanetary dust (IPD) in the optical and near-infrared, 
and thermal emission from the same IPD in the mid- and 
far-infrared. Historically, extensive zodiacal light obser- 
vations in the scattered sunlight regime were conducted 
from ground-based observations at high altitude sites 
in the 1960's and 1970's (Levasseur-Regourd & Dumot 
1980). However, these were restricted to optical wave- 
lengths because atmospheric OH airglow is much brighter 
than the zodiacal light in the near-infrared. Space-based 
platforms get above this contamination and provide pre- 
cision measurements of the diffuse sky brightness at near- 
infrared wavelengths (Murdock & Price 1985; Berriman 
et al. 1994; Matsuura et al. 1995; Matsumoto et al. 1996). 
Based on measurements from DIRBE on COBE^ phe- 
nomenological models for the IPD distribution have been 
constructed (Kelsall et al. 1998; Wright 1998) to estimate 
the zodiacal foreground for measurements of the extra- 
galactic background light (EBL; Wright (1998); Hauser 
et al. (1998b); Matsumoto et al. (2005)). 
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A continuous supply of IPD particles is necessary to 
replenish the zodical dust cloud, as they either fall into 
the Sun through Poynting-Robertson drag or leave the 
inner solar system due to radiation pressure. However, 
it still remains unclear whether the dominant source of 
IPD is asteroids or comets. Recent dynamical analyses 
suggest that IPD arises from a cometary origin (Liou et 
al. 1995; Ipatov et al. 2008; Nesvorny et al. 2010); dust 
trails of short-period comets show that the mass loss from 
cometary nuclei amounts to a third to a half of the re- 
quired supply of dust (Ishiguro & Ueno 2009). Addition- 
ally, particles captured in the stratosphere suggest that 
IPD from asteroids and comets occur with comparable 
abundance with cometary particles being slightly more 
common (Schramm et al. 1989). On the other hand, 
three dust bands in the IPD cloud were found by IRAS 
(Low et al. 1984) and their parent bodies were identified 
as being associated with certain asteroid families (Der- 
mott et al. 1984; Nesvorny et al. 2003). Age estimates 
of the last collisions in these families are consistent with 
such an association (Nesvorny et al. 2003). The discov- 
ery of these dust bands is strong evidence that there is 
an important asteroidal contribution to IPD. 

Because materials likely to be present in the IPD have 
distinct spectral features in the near-infrared, spectro- 
scopic measurement of the zodiacal light may help to 
determine the composition of IPD particles. Measure- 
ments based on bandpass photometry indicate that the 
color of the zodiacal light is redder than the solar spec- 
trum (Matsuura et al. 1995; Matsumoto et al. 1996), 
which implies that dust particles larger than 1 /im are 
mainly responsible for the zodiacal light (Matsuura et al. 
1995). Matsumoto et al. (1996) found a silicate-like fea- 
ture at around 1.6 jam from Infrared Telescope in Space 
(IRTS) data, and suggested that dust from S-type aster- 
oids are responsible for the zodiacal light. In addition, 
a silicate feature was found in mid-infrared spectrum of 
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Table 1 

Summary of the CIBER Observation Targets 
Target Time [sec] Altitude [km] R.A. Decl. X ^ I b Zenith Angle Solar Elongation 



Elat-10 


90 


- 109 


139 - 


Elat-30 


127 


- 142 


199 - 


Bootes-A 


162 


- 250 


245 - 


Bootes-B 


262 


- 358 


314 - 


NEP 


370 


- 426 


288 - 


SWIRE ELIAS-Nl 


440 


- 490 


215 - 



172 234.05 -8.32 233.76 10.71 

220 222.75 20.56 212.82 35.10 

310 218.63 34.84 201.23 46.28 

295 217.33 33.39 202.23 47.28 

232 270.29 65.88 - 90 

135 242.56 55.21 208.40 72.63 



357.23 36.63 44.4 102.0 

25.90 61.96 74.7 117.1 

58.59 66.68 84.5 119.3 
55.43 68.02 84.5 119.3 

95.60 29.69 42.7 89.9 
85.21 44.66 59.3 100.6 



Note. — The time column denotes the elapsed time after launch and (A,/3) and {l,b) denotes ecliptic and Galactic coordinates, respectively. 
All angles are given in degrees. 



the zodiacal light by IRTS (Ootsubo et al. 1998), In- 
frared Space Observatory {ISO; Reach et al. (2003)), and 
^ir^i?/ (Ootsubo et al. 2009). 

2. THE COSMIC INFRARED BACKGROUND EXPERIMENT 
2.1. Overview of CIBER 

The Cosmic Infrared Background ExpeRiment 
(CIBER) is a rocket-borne instrument to search for 
signatures from first-light galaxies in the near-infrared 
extragalactic background (Bock et al. 2006, Zemcov et 
al. 2010). CIBER consists of four optical instruments: 
two wide-field imagers (Sullivan et al. 2010); a Narrow 
Band Spectrometer (NBS; Renberger et al. 2010); and 
a Low Resolution Spectrometer (LRS; Tsumura et 
al. 2010). The LRS data from CIBER's first fiight 
are presented in this paper. These instruments are 
mounted to a common optical bench which is cooled 
by liquid nitrogen to reduce thermal instrumental 
emission. A Terrier-Black Brant IX rocket carrying the 
CIBER payload was successfully launched from White 
Sands Missile Range on 2009 February 25, achieving an 
apogee of ~330 km, and providing more than 420 s of 
astronomical data. After the fiight, the instrument was 
successfully recovered for future fiight s. 

The choice of our science fields is shown in Figure 1 
and Table 1. Since the Bootes and SWIRE ELIAS-Nl 
fields were well studied with Spitzer (Eisenhard et al. 
2004; Lonsdale et al. 2003) and the north ecliptic pole 
(NEP) field is weh studied with AKARI (Matsuhara et 
al. 2006), these fields were chosen for the CIBER imaging 
instruments. The Elat-10 and Elat-30 fields were chosen 
specifically to measure the dependence of the zodiacal 
light intensity on ecliptic latitude, while the remaining 
science fields were chosen at higher ecliptic latitudes to 
reduce the brightness of the zodiacal foreground for the 
EBL science. However, as the zodical light intensity is 
never small, the EBL fields can also be used to constrain 
the zodical light. Each of these fields was observed in a 
series of 25 seconds reset intervals. As the Bootes field 
ancillary coverage is very wide, this field was observed 
in two pointings, A and B, with each observed for about 
90 seconds. During the last 50 seconds on Bootes-B, 
the cold shutters on the both spectrometers were closed 
to obtain dark frames for absolute spectroscopy. The 
pointing stability in 25 seconds was less than 3 arcsec. 

2.2. The Low Resolution Spectrometer 

The LRS is designed to obtain absolute spectra of 
the astrophysical sky brightness in the near-infrared. 
Although primarily designed to measure the absolute 



350 




Time [sec] 



Figure 1. Instrument altitude and science observation field se- 
quence as a function of time for CIBER's first fiight. Reference 
calibration lamps and cold shutters in front of the focal planes 
were operated in fiight to allow calibration transfer and absolute 
dark level reference. 

sky brightness from extragalactic sources, including first- 
light sources present during re-ionization, the LRS also 
measures the zodiacal light foreground. 

The specifications of the LRS are summarized in Ta- 
ble 2. The LRS is designed to cover the wavelength range 
0.75-2.1 fim with a spectral resolution A/ A A varying be- 
tween 15 and 30. The LRS employs a 5 cm refractive tele- 
scope and has a cold shutter just before the focal plane 
detector to monitor the detector dark current for abso- 
lute spectroscopy. The shutter is closed prior to launch, 
on ascent, during fiight, and on descent. The LRS optics 
form an initial focus at a field stop, where a mask with 
five equally spaced slits is placed. The optics then re- 
lay the field stop back to a parallel beam where a prism 
disperses incident light perpendicular to the slit mask 
direction. Finally, this parallel beam is refocused on a 
256x256 substrate-removed HgCdTe PICNIC detector 
array; the slit mask is imaged to five 2.8 arcmin x 5.5 
degree strips on the sky at the focal plane. A calibra- 
tion lamp is operated both before and after astronomical 
observations to confirm the stability of the instrument 
responsivity during the fiight: the detector was stable to 
better than 4% during the fiight. This calibration lamp 
also serves as a transfer standard between the fiight data 
and laboratory-based calibration measurements. 

The wavelength and responsivity calibration of the 
LRS were measured in the laboratory before fiight. For 
the wavelength calibration, radiation from a halogen 
lamp was dispersed by a monochromator, the output of 
which is fiber coupled to an integrating sphere produc- 
ing an approximately solar spectrum, aperture-filling in- 
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Table 2 


Specifications of the Low Resolution Spectrometer 


Characteristic 


Value 


Optics 


14 lenses, 2 filters, 1 prism, 5 slits 


Aperture 


50 mm 


F number 


2 


FOV 


5.5 degrees along a slit 


Pixel size 


1.4 arcminxl.4 arcmin 


Slit 


2 pixels X 236 pixels 


Wavelength range 


0.75-2.1 /xm 


Wavelength resolution 


A/AA=15-30 


Optics efficiency 


0.8 


Detector 


256x256 substrate removed PICNIC 


Detector QE 


0.9 


Dark current 


0.6 e-/s 


Readout noise 


26 e- CDS 



put source. The LRS observes the integrating sphere 
while the wavelength is stepped from 0.75 to 2.1 /im 
in 0.005 fim increments. Additionally, the responsivity 
of the LRS was calibrated using the halogen lamp and 
an optical/near-infrared tunable laser covering the range 
from 0.7 to 1.6 /im. Each of these sources is coupled to 
an integrating sphere and is measured simultaneously by 
both the LRS and absolutely calibrated detectors (SIR- 
CUS; Brown et al. 2004). These measurements are used 
to accurately calibrate both the absolute responsivity 
and the channel-to-channel relative responsivity of the 
LRS. 

Further details of the instrument including the calibra- 
tion and flight performance are summarized in Tsumura 
et al. (2010). 

3. DATA REDUCTION 

3.1. Obtained Images 

Figure 2 shows an LRS image obtained during 
CIBER's first flight. The five vertical slits are dispersed 
horizontally so that the wavelength increases from left to 
right in this image. The horizontal lines are caused by 
stars which happened to illuminate the slit, while the five 
strong vertical lines are caused by stray thermal emission 
from the hot rocket skin prominent at long wavelengths 
(see Section 3.2). The central slit has a small notch lo- 
cated at the bottom center used to facilitate laboratory 
testing. A masked region exists at one edge of the ar- 
ray which can be used to monitor the presence of diffuse 
stray light. 

3.2. Thermal Stray Light 

The five strips visible at longer wavelengths in Figure 2 
are caused by bright thermal emission from the sound- 
ing rocket skin which was heated by air friction during 
the powered ascent. This thermal emission couples to 
the LRS optics via diffuse scattering from the instrument 
baffle placed before the first aperture. By fitting a black- 
body curve to the thermal spectrum in the LRS images 
taken during the time when the shutter door was closed, 
the temperature of the rocket skin is estimated to have 
peaked at ~400 K in flight. Spillover from this stray 
thermal emission from the adjacent slit contaminated 
the short wavelengths regions (<0.82 /im) of the four 
rightmost slits, so only data from the leftmost slit are 
used to derive the four spectral data points shorter than 
0.82 /im. Contamination by longer wavelength spillover 



Figure 2. Image of the NEP field obtained during 25 seconds of 
integration of the LRS during CIBER's first fiight. The dispersion 
direction is horizontal and the wavelength increases from left to 
right in the images, while the large intensity in the five perpendic- 
ular strips is caused by scattered thermal emission from the hot 
rocket skin. The central slit has a small notch located at the bot- 
tom center to facilitate laboratory testing; the edge of the array 
is not illuminated by the optics and can be used as a monitor of 
diffuse stray light falling on the detector. The left hand panel in- 
cludes the raw spectra of stars which happened to fall on the slits 
as horizontal bands of emission, and the detected stars are masked 
as shown at the right hand panel. The averaged spectrum of the 
diffuse sky is obtained by combining these source-masked spectra. 

within each slit region was measured after the fiight to 
be less than 0.4% for A <1.4/im (Tsumura et al. 2010), 
which is negligibly small for the science described in this 
paper. In addition, we investigated the masked region 
of the astronomical data and found a low photo current 
level (<2 e~/sec corresponding to <100 nW/m^/sr), a 
small residual which we subtracted. There is no evidence 
that stray thermal emission affects the results presented 
here at shorter wavelengths, and we simply omitted data 
at A >1.8 /im. 

3.3. Airglow Emission 

A second emissive component not associated with the 
astronomical sky is observed to have exponential time 
and altitude dependence as shown in Figure 3. Based 
on its spectrum, which has prominent emission features 
at 1.1 and 1.6 /im, we attribute this component mainly 
to OH molecules from either exospheric atmospheric air- 
glow or dissociated water vapor which outgasses from 
the pay load early in the flight. The spectral shape of 
the airglow emission is observed to be the same during 
ascent and descent, but as shown in Figure 3 the abso- 
lute brightness of the airglow emission on ascent is much 
brighter; this suggests an enhancement by dissociated 
water outgassing from the pay load early in the flight. 
This airglow emission is bright in the Elat-10 and Elat- 
30 field data observed early in the flight, and is essen- 
tially negligible for the remaining fields; Figure 3 shows 
the brightness of the airglow emission at 1.61 /im in each 
field derived by subtracting the NEP brightness. 

As the Elat-10 field was observed earliest in the fiight 
and at the lowest altitude where we would expect this 
airglow emission to be brightest, these data are used to 
derive a template spectrum of the contaminant. In order 
to subtract the effect of the airglow emission, a spectral 
model was derived. The total signal / can be modeled as 

/ = A(A)+B(A)exp(-t/r), (1) 

where A{X) is the astrophysical sky brightness, B{X) is 
the airglow emission, and r is the timescale of the decay. 
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Figure 3. (a) Brightness of the airglow emission as a function 
of instrument altitude at 1.61 jim where the airglow emission is 
brightest. The Elat-10, Elat-30, and SWIRE data are separated 
into several subsets to allow differentiation of the time and altitude 
dependence of the airglow signal. The airglow emission appears to 
be negligible at altitudes greater than 250km. (b) Brightness of 
the airglow emission at 1.61 /im which scales exponentially with 
time during the ascent phase of the flight. A similar exponential 
variation is observed on descent but with a smaller amplitude. The 
decay timescale r was 14.2 seconds from only the Elat-10 data, 
while we obtained r=19.9 seconds for the entire data sequence. 
We used r=14.2 seconds for deriving the absolute brightness of 
the airglow emission at Elat-10 fleld. 

Here we employ a fit to airglow in time not in rocket's 
altitude because the airglow emission on ascent is domi- 
nated by the outgassing, and the outgassing should follow 
the exponential decay in time rather than in altitude. To 
increase the constraining power of the data, the Elat-10 
data is reduced into four and Elat-30 into three sequen- 
tial spectra. Because the sky the LRS observes during 
the observation of a single field is static, any change in the 
mean absolute brightness must be caused by time varia- 
tion in the airglow component. The spectral shape of the 
airglow component B{X) is derived by a difference of the 
brightness between the first block and the last block in 
Elat-10. The decay timescale r was also derived with a 
fit to the data in Figure 3 and we obtained r = 14.2 sec- 
onds from Elat-10 data. Combining the differencing and 
the obtained decay timescale, we derived the absolute 
brightness of the airglow emission at Elat-10. 

The absolute brightness of the airglow emission at 
Elat-10 derived using this method was checked by 
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Figure 4. Subtraction of the airglow emission for the Elat-10 and 
Elat-30 fields observed early in the ascent phase where airglow is 
brightest. The blue points are "original" spectra of the diffuse sky 
from about 20 seconds integration, which include both the airglow 
emission and astrophysical sky brightness. The green points are 
spectra of the airglow emission at each field derived by the method 
described in Section 3.3. The red points show the residual spectra 
obtained after subtracting the airglow emission. 

comparing the color of the residual spectrum between 
1.61 /im (maximum airglow emission) and 1.29 /im (lo- 
cal minimum airglow emission) to the spectrum of the 
other fields. We confirmed that the color of the resid- 
ual spectrum of Elat-10 was consistent with the color of 
the NEP and Bootes spectra to within 5%. The color of 
the SWIRE ELIAS N-1 field was 10% higher than the 
others because SWIRE was observed at the end of the 
flight sequence and the airglow emission began to rise 
again on descent. Therefore, we attribute a 5% error to 
our estimation of the airglow emission based on the color 
of the residuals. The absolute brightness of the airglow 
emission at the Elat-30 field was estimated by the same 
color estimation method and subtracted, that is, the air- 
glow intensity at Elat-30 was set to make the color of the 
residual spectrum between 1.61 /im and 1.29 /im same as 
NEP and Bootes. Figure 4 shows the pre- and post- 
subtraction spectra with the airglow emission. We are 
confident that our main scientific result on the detection 
of the 0.9 /im absorption feature is not strongly affected 
by airglow emission. 
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Table 3 

Comparison of Total and Integrated Galactic Star Light at J-band 



Sky Brightness 


Elat-10 


Elat-30 


Bootes 


NEP 


SWIRE 


(a) CIBER/LRS intensity at 1.21 /im 


792.09 


465.29 


406.42 


417.74 


388.17 


(b) Integrated starlight (J-band >13 mag) 


46.43 


20.83 


18.06 


44.71 


27.31 


(a) / (b) [%] 


5.86 


4.48 


4.44 


10.70 


7.03 



Note. — Values in (a) and (b) are shown in nW/m^/sr. 



3.4. Point Source Detections 

The flight data include spectra of stars which fell on 
the slits as bright horizontal lines as shown in Figure 2. 
All detected horizontal lines were collected and masked 
when deriving the diffuse sky spectrum. In this proce- 
dure, the images of point sources with mvega(-^-band)>13 
(Fo(Vega) = 1630 Jy) were masked. Remaining sources 
are also removed when we average the source-masked 
spectra with 3a clipping to derive the diffuse sky spec- 
tra. To estimate the contribution to the total sky bright- 
ness from Galactic stars fainter than the LRS detection 
limit, all J-band stars in the 2MASS point source cata- 
log (Skrutskie et al. 2006) from 2° x 2° regions centered 
at the each fleld are summed; the lOa detection limit in 
2MASS is 15.8 mag (Vega magnitude) in J-band. To un- 
derstand the completeness of the 2MASS-determined in- 
tegrated starlight, the value determined from the 2MASS 
collection is compared with that computed using the 
SKY model which is designed to provide predictions for 
infrared Galactic star counts (Cohen 1994), and it is 
found that the calculated 2MASS integrated intensity 
was ~90% of the SKY model estimation. The SKY 
model takes into account the flux from stars down to 
a much deeper limit of ~30 mag which accounts for this 
difference. We therefore corrected the calculated 2MASS 
integrated star light based on the SKY model estimation 
to account for stars fainter than 15.8 mag, and the re- 
sults are summarized in Table 3. The contribution to the 
sky brightness from residual Galactic stars was ~ 5% be- 
sides NEP, where the stellar contamination was largest 
(~ 10%) because of its low Galactic latitude. 

The spectra of bright reference stars are used to check 
both the relative in-flight channel-to-channel response 
of LRS and the ground-based wavelength calibration. 
Because the absolute pointing of the instrument is not 
known with sufficient precision, it is not possible to de- 
termine exactly where on the slits bright stars fell and 
so we are not able to use this information for absolute 
brightness calibration. Figure 5 shows the measured LRS 
spectrum of 42 Draconis (42Dra, a magi:f=2.2, K1.5III 
star), the brightest star observed with the LRS during 
the flight. Also shown for comparison are the spectrum 
of 42Dra from the Pulkovo Spectrophotometric Catalog 
(Alekseeva et al. 1996), the 2MASS fluxes (Skrutskie et 
al. 2006), and a model of a K1.5III star a-Boo (Cohen et 
al. 2003). The LRS spectrum of 42Dra is 2.1 times fainter 
than the absolute brightness of the Pulkovo value; this 
difference is due to partial obscuration of the star by the 
slit mask. Although the model spectrum does not accu- 
rately reproduce spectral features observed with the LRS 
because the accuracy of small spectral features in the 
model spectrum is not very good, the LRS and Pulkovo 
spectra both exhibit such features in their regions of over- 
lap. Accounting for the accuracy of the model, the scaled 



LRS spectrum is consistent with the expected spectrum, 
conflrming the relative response calibration obtained in 
the laboratory. 

2.5 X 1 0"^ I I 
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Figure 5. Measured LRS spectrum of 42Dra; also plotted are the 
equivalent spectrum from the Pulkovo Spectrophotometric Cata- 
log (Alekseeva et al. 1996), a model spectrum (Cohen et al. 2003), 
and photometric data from 2MASS (Skrutskie et al. 2006). Be- 
cause 42Dra was partially obscured by the slit mask, the curve 
has been scaled by a factor of 2.1 to normalize it about the ab- 
solutely calibrated spectra. The LRS spectrum agrees with the 
other data, confirming the relative channel-to-channel calibration 
obtained in the laboratory. The periodic spectral features evident 
in the Pulkova and LRS spectra are not reproduced by the model 
because of its less accuracy for the small spectral features. 

4. RESULTS 
4.1. Sky Brightness 

Figure 6 shows the spectra of the sky brightness for 
fields observed during CIBER's first fiight. The Elat- 
10 and Elat-30 plots show the spectra after subtract- 
ing airglow emission; airglow emission in the other fields 
is negligible. The LRS sky spectra show the total 
brightness from zodiacal light, Galactic stars, and the 
EBL, and are consistent with previous observations with 
COBE/DIRBE (Hauser et al. 1998b) and IRTS (Mat- 
sumoto et al. 1996) in the region of overlap. The spectral 
shapes of the sky brightness in Figure 6 are very similar 
to each other, irrespective of ecliptic latitude. A previ- 
ously unreported broad absorption feature centered at 
approximately 0.9 fim is present in each of these spectra. 

Figure 7 shows the ecliptic latitude dependence of the 
zodiacal light calculated using the DIRBE all-sky zodia- 
cal dust model (Kelsall et al. 1998) and the observed sky 
brightness after removing the integrated Galactic star 
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light at 1.25 /im. Except at high echptic latitude, the 
observed LRS brightness agrees with the ecliptic latitude 
dependence estimated by the model, which shows again 
that our observation is consistent with the previous ob- 
servations. We note that the zodical light model departs 
from the data at higher ecliptic latitude. 
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Figure 6. Measured spectra of the near-infrared sky brightness. 
Colored circles indicate the CIBER/LRS data used in this study, 
open squares indicate IRTS data averaged at 72° -73° ecliptic lat- 
itude from Matsumoto et al. (1996), and purple squares indicate 
darkest DIRBE data from Hauser et al. (1998b). The error bars 
denote the combination of statistical and systematic error due to 
subtracting airglow emission. The absolute calibration error is es- 
timated as ±10% for LRS (Tsumura et al. 2010), ±5% for IRTS 
(Noda et al. 1996), and ±1.6% for DIRBE (Hauser et al. 1998b). 
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Figure 7. Ecliptic latitude dependence of the DIRBE all-sky zo- 
diacal dust model (Kelsall et al. 1998) (open squares) and the LRS 
sky brightness (circles) after removal of the integrated Galactic 
starlight at ~ 1.25 /im. The impact from the difference of the solar 
elongation or ecliptic longitude should be small, because the range 
of them in our observation fields is narrow (<30°; see Table 1) to 
see the ecliptic latitude dependence of the zodiacal light. 



4.2. Zodiacal Light Spectrum 

To separate residual astrophysical components from 
zodiacal light, the sky spectra observed in two fields 
are differenced; as shown in Figure 8, various combi- 
nations of such field differences show a similar spec- 
tral shape suggesting that the zodiacal light spectrum 
is largely isotropic. A joint CIBER-IRTS zodiacal light 
spectrum is produced by normalizing both data sets to 
the DIRBE all-sky zodiacal dust model (Kelsall et al. 
1998) at 1.25 /im and 2.2 /im using the bandpass response 
of DIRBE filters (Hauser et al. 1998a). Figure 9 shows 
the derived joint CIBER-IRTS zodiacal light spectrum 
with the solar spectrum (Gueymard et al. 2002). The 
measured zodiacal light spectrum has a redder color than 
the solar spectrum, and is generally similar to the unsub- 
tracted spectra shown in Figure 6, though with much less 
signal to noise due to the differencing. Hereafter, LRS 
data at >1.4/im were removed because their error bars 
were too large because of the subtraction of the large 
airglo" 
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Figure 8. Zodiacal spectra from various field difference combi- 
nations. The figure (a) shows the difference between Elat-10 and 
the science fields, while the figure (b) shows the difference between 
Elat-30 and the science fields. The curves are scaled to separate 
them for visualization. The solid lines denote the observed solar 
spectrum (Gueymard et al. 2002) normalized at 1.08 /im. The dif- 
ference spectra are noisier than the raw spectra as expected, but 
remove EBL and partially Galactic emission and show the absorp- 
tion feature, confirming it is of zodiacal origin. 
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Figure 9. Zodiacal light spectrum. The filled circles show the 
CIBER/LRS data obtained by taking an average of the differences 
between Elat-10 and the other fields. The IRTS data (Matsumoto 
et al. 1996) are plotted in open squares and the curve denotes the 
observed solar spectrum (Gueymard et al. 2002) normalized to the 
IRTS data at 1.83 /im. The GIBER and IRTS data are normalized 
to the brightness expected at the Elat-30 field (see Table 1 for 
specific coordinates) using the DIRBE zodiacal light model plotted 
in open diamonds (Kelsall et al. 1998). Error bars on the LRS data 
denote the combination of statistical error and systematic error in 
subtracting airglow emission. Error bars on the DIRBE data show 
the CIBER-IRTS relative calibration error effecting the relative 
color calibration between GIBER and IRTS which is estimated to 
be 3.5% from the color difference of the two different zodiacal light 
models (Kelsall et al. 1998; Wright 1998). The relative scaling 
error due to the zodiacal light model is 10% estimated from the 
difference of the absolute zodiacal intensity in the two models. 



4.3. Reflectance of the Interplanetary Dust 

Figure 10 shows the IPD reflectance spectrum, de- 
rived by dividing the LRS and IRTS spectra by the solar 
spectrum. Note that the reflectance of zodiacal light as 
viewed from Earth is dominated by the IPD near Earth's 
orbit because the radial density profile of IPD results 
in a contribution to observed intensity that falls off as 



-2.5 



(Leinert et al. 1998). The contribution of dust 
scattering at distances greater than 2.5 AU to the in- 
tegrated zodiacal light intensity is estimated to be less 
than 3% based on this radial density profile. The re- 
flectance of an S-type asteroid 25143-Itokawa (Binzel et 
al. 2001), a C-type asteroid 1-Ceres (Bus & Binzel 2002), 
and a comet 9P/Tempel-1 (Hodapp et al. 2007) are also 
shown in Figure 10 for comparison. The C- and S-type 
asteroids are respectively the first and second most com- 
mon types of asteroid in the asteroid belt. Unfortunately, 
due to the difficulty of such observations no spectra of 
old cometary dust trails at near-infrared wavelengths are 
available. Instead we show the spectrum of eject a ma- 
terial from comet 9P /Temple- 1 after the impact of the 
deep impact probe (Hodapp et al. 2007); this spectrum 
is of "fresh" dust and other materials from the cometary 
nucleus and so is not an ideal comparison to the zodiacal 
reflectance spectrum. 
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Figure 10. Reflectance of IPD normalized at 1.83 /im. LRS and 
IRTS results are plotted together with example spectra of comets 
and asteroids from the literature. The reflectance spectra of the 
LRS (fllled circles) and IRTS (open squares) were obtained by di- 
viding the zodiacal light spectra by the solar spectrum. Curves 
show the reflectance of an S-type asteroid 25143-Itokawa from 
Binzel et al. (2001), a C-type asteroid 1-Ceres from Bus & Binzel 
(2002), and a comet 9P/Tempel-1 from Hodapp et al. (2007), re- 
spectively. 



5. DISCUSSION 

5.1. Composition of the Interplanetary Dust 

We first detected the broad absorption feature in the 
IPD reflectance spectrum at 0.9 /im in Figure 10. This 
absorption feature is common in certain types of as- 
teroids and meteorites, and it is thought to be caused 
by silicate compounds such as pyroxene and/or olivine, 
which are the principal mineralogical constituent of IPD 
(Bradley 2003). Variation of this feature over a wide 
range of ecliptic latitude (10-90°) is not detected, which 
is consistent with mid-infrared observations (Reach et al. 
2003). 

Ordinary chondrites have known absorption bands at 
1 /im and 2 /im, evidence for the presence of pyroxene 
and olivine. The spectrum of 25143-Itokawa has been 
described as being similar to that of ordinary chondrite 
meteorites with additional reddening due to space weath- 
ering (Binzel et al. 2001), a process that alters the state 
of surface materials exposed to a space environment. 
The exposed surface develops radiation damage, vapor- 
and sputter-deposited coatings, and melt products in- 
cluding agglutinates. The reddening effect particular 
to space weathering is caused by a vapor coating con- 
taining nanophase-reduced iron (npFe^) particles around 
each dust particle (Pieters et al. 2000; Nimura et al. 
2008). The spectrum of 25143-Itokawa is well matched 
to a model including 0.05% npFe^ reddening of ordinary 
chondrite (Binzel et al. 2001), a value common to S-type 
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asteroids (Binzel et al. 2001; Hapke 2000). The simi- 
larity of the shape of the IPD reflectance measured by 
LRS to that of S-type asteroids is evidence that this pro- 
cess is relevant to the IPD responsible for producing the 
zodiacal light in the near- infrared. 

Matsumoto et al. (1996) reported some similarity be- 
tween the spectrum of zodiacal light and S-type asteroids 
from the IRTS data. Our result is also consistent with the 
finding of a silicate feature in the mid-infrared zodiacal 
light spectrum (Ootsubo et al. 1998, 2009; Reach et al. 
2003), though the peak position of the absorption band 
in the zodiacal light spectrum is slightly shifted from that 
of S-type asteroids. This difference may be explained by 
differences in the chemical composition and/or some de- 
pendence on the size of the IPD particles. 

5.2. Distribution of the Interplanetary Dust 

As it is difficult to formulate a mechanism which pro- 
pels asteroidal dust out of the ecliptic plane, recent dy- 
namical analysis examining the ecliptic latitude depen- 
dence of the zodiacal emission at mid-infrared wave- 
lengths suggests that the IPD arises mainly from a 
cometary population (Nesvorny et al. 2010). Though the 
zodiacal emission is biased toward low-albedo dust from 
comets and C-type asteroids with albedo ~0.05 (Tedesco 
et al. 1989), the scattered light seen by the LRS may 
be biased by high- albedo S-type dust with albedo ^0.2 
(Tedesco et al. 1989). The result presented here can be 
consistent with the resent result favoring a cometary ori- 
gin if high-albedo asteroidal dust is concentrated near 
the Earth. Figure 11 shows the histogram of asteroids in 
the asteroidal belt as a function of their inclination, in- 
dicating that asteroids are distributed over 10°-15° incli- 
nation at half maximum. Assuming that the asteroidal 
dust generated in the asteroidal belt falls into the sun 
purely by the Poynting- Robertson effect, the estimated 
scale height of the dust cloud from the asteroid belt at 
1 AU is ~0.2 AU. This thickness is comparable with the 
estimated IPD cloud scale height of the DIRBE all-sky 
zodiacal dust model (Kelsall et al. 1998). Cometary dust 
may be distributed at distances greater than 0.2 AU from 
the Earth, if the ecliptic latitude dependence of the zodi- 
acal light is masked by the nearest high-albedo asteroidal 
dust. 

The spatial distribution of asteroid types within the as- 
teroid belt is non-uniform. Although S-type asteroids are 
the second largest population within the asteroidal belt, 
they predominantly occupy the inner portion of the belt. 
The relative fraction of S-type asteroids at heliocentric 
distances out to 2.5 AU is no less than 60% (Mothe-Diniz 
et al. 2003). Using this fractional abundance and assum- 
ing that dust production rate of each asteroid type is pro- 
portional to the square of the number density, the dust 
ratio of S-type to C-type asteroids is expected to be 5:3. 
Furthermore, the C-type asteroids have a smaller albedo 
(~0.05) compared to S-type asteroids (~0.2) (Tedesco et 
al. 1989). The combination of asteroid type ratio and the 
albedo ratio suggests that S-type dust dominates over C- 
type dust in the observed zodiacal light spectrum by a 
fraction of 6 and 7. 

These estimates are qualitatively consistent with the 
interpretation that the CIBER/LRS zodiacal light spec- 
trum shows the existence of the S-type asteroidal dust 
particles, and that they are largely responsible for the 



scattered zodiacal light spectrum observed from Earth. 
Our result, however, does not rule out an important con- 
tribution from comets. We are unable to state the exact 
dust fraction from comets versus asteroids as the high 
albedo dust near the Earth's orbit dominates the ob- 
served spectrum. Furthermore there is still the possi- 
bility that cometary dust might also show the pyrox- 
ene/olivine feature at 0.9 /im. Pyroxene and olivine 
have been confirmed in the mid-infrared spectrum of the 
cometary dust (Lisse et al. 2006) and in the cometary 
dust samples returned by the Stardust spacecraft (Zolen- 
sky et al. 2006), but there is no evidence of these ab- 
sorption features in the near-infrared spectrum. Spec- 
troscopic observations of the cometary dust trails near 
1 /im will be necessary to constrain the cometary contri- 
bution to the observed LRS spectrum, and to estimate 
the fraction of dust from comets in the IPD. A measure- 
ment of near-infrared and mid-infrared spectral features 
obtained in a deep space mission traveling to the outer so- 
lar system could map the composition of IPD throughout 
the solar system, and trace its origins to asteroidal and 
cometary parent populations (Matsuura 2002; Cooray et 
al. 2009). 
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Figure 11. Histogram of the asteroids in the asteroidal belt 
as a function of inclination of the orbit. The complete 
data set of the numbered asteroids from the JPL Web Site 
(http://ssd.jpl.nasa.gov/7sb_elem) was used. 
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